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ABSTRACT 

A new equivalence relation, named relation of 'similarity' is defined and ap- 
plied in the restricted three-body problem. Using this relation, a new class of 
trajectories (named 'similar' trajectories) are obtained; they have the theoretical 
role to give us new details in the restricted three-body problem. The 'simi- 
lar' coordinate systems allow us in addition to obtain a unitary and an elegant 
demonstration of some analytical relations in the Roche geometry. As an exam- 
ple, some analytical relations published in Astrophysical Journal by Seidov in 
2004 are demonstrated. 

Subject headings: Restricted problems: restricted problem of three bodies — Stellar 
systems: binary stars 
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1. Introduction 

In the frame of the restricted three-body problem, the Roche geometry is a fundamental 
notion, and it is much studied, using for this aim different coordinate systems. 

The usual way of treating the problem of motion of test particle and the Roche 
geometry in the gravitational field of a binary system consists in introducing a coordinate 
system (xyz), rotating jointly with the components. The x-axis passes through the centers 
of both components, and the y-axis is situated in the orbital plane; but there are more 



possibilities to locate the origin of the coordinate 



located in the mass center of the binary system ( Moulton 



more m assive star (IKopal 



papers ( iHuang 



1967 



1978 



1989 



Kruszewski 



Rov 



1988 



^ system. '. 


-^or ex< 


(Moulton 


1923). 


Szebehelv 


1967) 



or ex ample the origin can be 



923), or in the center of the 



19671 ). But there are some 



19631 ). where the location of the origin of the coordinate 



system is not precisely indicated. For example Huang (1967) wrote: "Thus, if we denote 
1 — fi as the mass of one component, fi will be the mass of the other. Let us further choose 
a rotating (x,y,z) system such that the origin is at the center of the 1 — /i component, the 
X-axis points always towards the /x component, and the xy plane coincides with the orbital 
plane." 

This way to locate the origin of the coordinate system is not an ambiguous one, but it 
offers the opportunity of the question: What happens with the equations of motion of the 
test particle and with the geometry of the equipotential curves, in the restricted three-body 
problem, if the origin of the coordinate system is taken in the center of the less massive 
star. 

The aim of this paper is to answer the question above, by introducing a new notion: 
the 'similar' coordinate systems. 



In order to do so, a binary relation is created, which is denoted by the author as 
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relation of 'similarity'. Then, using this relation, the 'similar' coordinate systems, and 
the necessary 'similar' parameters and physical quantities are defined, obtaining 'similar' 
equations of motion, and 'similar' equipotential curves. 

The conclusion of this study is that the use of 'similar' coordinate systems in the 
restricted three-body problem allows us to complete the traditional study of the Roche 
geometry with some new features. 

2. Relation of 'similcirity' 

In this article we shall write 'similar' (not similar), because we intend to use this word 
as the name of a new mathematical relation, and not as an adjective. 

Definition: Two or more mathematical objects are in the 'similarity' relation in 
connection with a given definition (V), if the objects are completely characterized by V. 

For example: 

1. in algebra: The numbers xi = 1 and X2 = —3 are 'similar' in connection with the 
definition: x is the solution of equation: x^ + 2x — 3 = 0, x G Z. 

2. in geometry: If we consider the definition: P is the point situated into a given plane 
(tt), at the distance r = 1 from the given point A, A & (tt), then all the points of the 
circle having the center A and radius r — 1 are 'similar'. 

3. in astrophysics: In the frame of the circular restricted three-body problem we define 
a comoving coordinate system situated in the orbital plane, having the origin in the 
center of one component of the binary system and the abscissa 's axis pointing to the 
other component. Then the coordinate systems xSiy and x'S2y' are 'similar' (see 
Figure 1). (In many books or articles of astronomy, this is the coordinate system's 
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definition used to study the restricted three-body problem (iHuang 1119671 ) . but usualy 
only xSiU coordinate system is considered. In 1963, Kruszewski wrote: "The center 
of coordinates is p laced at the center of the (arbitrarily chosen) primary component" 



( iKruszewski 



19631 ). This means that if the origin of the coordinate system is taken 
into the center of the secondary component, the results will be similar. In this article 
we try to find what this similarity implies.) 



The study of the restricted three-body problem in 'similar' coordinate systems follows 
the classical algorithm, but some typical peculiarities appear. So, the use of the 'similar' 
coordinate systems impose the use of some physical and geometrical 'similar' quantities: 

- 'similar' mass ratios q and q' 

- 'similar' distances ri, r2, r'^, and r'2 

- 'similar' initial velocities fox, voy, v'q^, and Vq^ 

- 'similar' initial positions Xq, yo, x'q, and z/'q. 

Of course, as in algebra the 'similar' solutions of a polynomial equation are connected by 
the relations of Viete, the two 'similar' coordinate systems are connected by the equations 
of coordinate transformations (see section 4). 

It is easy to verify that the relation of 'similarity' is reflexive, symmetric, and 
transitive. That means that the relation of 'similarity' is an equivalence relation. 

Remark: 



1. The role of the definition into the establishment of a 'similarity' is huge. So, if we 
consider the definition: x is the solution of equation: -|- 2a; — 3 = 0, a; G N, the 
numbers Xi = 1 and X2 = —3 are not 'similar'. 
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2. The solutions of a problem (including the repeated ones) are in relation of 'similarity', 
because these can be defined as satisfying the same problem. 



3. 'Similar' coordinate systems, 'similar' parameters and physical quantities 



19671 ). we will 



In the frame of the circular, restricted three-body problem (iSzebehelyl 
consider Si and 5*2 the components of a binary system, mi and m2 their masses and Li, 
i = 1,5 the Lagrangian points. Due to the normalization (see section 4), the distance 
between and 5*2 is equal to 1. We will consider a rectangular coordinate system, so that 
one component of the binary star has coordinates (0,0,0), the other one has coordinates 
(1,0,0), and the angular Keplerian velocity Uk has components {0,0, Uk)- We observe 
immediately that there are two coordinate systems: (xSiyz) and {x'S2y'z'), which can be 



built. These are ' similar ' coor dinate systems (see Figure 1). As it is well-known (Kopal 



19781 ) p. 327-328, (IKopall Il989h p. 15-16 the mass ratio q is the main parameter which 
describes the Roche geometry. If we denote q = (mass which isn't into the origin) / (mass 
which is into the origin), we have q = — and q' = — which are 'similar' parameters. 



We denote r\ = distance of the infinitesimal mass (jSzebehely 



m2 



19671 ) to the origin of the 



coordinate system, and r2 = distance of the infinitesimal mass to the star which is not into 




Fig. 1. — The 'similar' coordinate systems with origin in 5*1 and 5*2 
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the origin. Therefore ri with r[ and r2 with r2 are 'similar' distances (see Figure 1). 

We have "'similar' velocities v = ^ and v' = and 'similar' accelerations a = 
and a' = 



4. 'Similar' equations of motion 



The forces w hich act on infinitesimal mass are Fg^avi, Fgrav2, Fcentrif, and Fcor 



iolis 



(jSzebehelyl 119671 ) p. 590). Their expressions in (xSiyz) and {x'S2y'z') coordinate systems 
are: 



gravl 



J 



k 



grav2 



Gm2m3[x — {Ri + -R2)] ^ Gm2m3y -t Gm2m^z 



k 



(2) 



^centrif — 1^3 



[x - Ri)i + yj 



(3) 



.dy^ dx^ 

^Coriolis — ^f^Z 1 — -^3 



and respectively: 



(4) 



grav\ 



Gm2'm^x' Gm2Tn^y' -t Gm2m^z' 



(5) 



grav2 



GmiTJi^lx' — [Ri + R2)] ^ Girtim^y' -r Gmim^z' 



k 



(6) 



_ 2 

-f^ centrif — ^k 



[x' - i?2)^''+ y'l 



(7) 
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F'corioiis = -2m3 Uk { —i - —j 



(8) 



where G is the gravitational constant. 



Fcorioiis and F'corioUs have opposite signs because of the orientation of vectors cj^ x w 
and u'k X v' (see Figure 2). 

By consequence in the (xSiyz) coordinate system the vectorial equation of motion is: 

^3 Fgravl ~\~ Fgrav2 ~\~ Fcentrif ~\~ Fccyriolis (9) 

and in the {x'S2y'z') coordinate system: 



Tfl^ CL F gravl ~l~ F grav2 ~l~ F centrif ~l~ F Coriolis 



(10) 



We shall use a special unit system (IKopall 119781 ) p. 318: we choose for the mass unit 
the sum of the masses of the components of the binary system, for the distance unit the 
distance between the centers of the components, and for the time unit the reciprocal of the 
angular Keplerian velocity. In that case, the orbital period will be P = 27r, and G = 1. 



Then mi = = Ri 



' R2 



1+g 



1 











/ 






7\ 

/ 1 


/ 

/ 






/ 1 

/ 








1 U 








/ 

/ 

/ 

/ 




Fig. 2. — Explanation for the opposite signs of Coriolis force in the two 'similar' coordinate 
systems . 
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Using the same unit system, the scalar equations of motion in the (xSiyz) coordinate 
system become: 

f^_2^^x - ^ - 1) (ID 

df^ dt 1 + q {l + q)rf {1 + q)r^ 



^ + 2—= ^ iy (-,0) 

dt^ dt ^ (1 + q)rf (1 + q)ri ^ ' 



d z z q z 



d'^z' z' q' z' 



dt^ (1 + (1 + r/)r;3 ' 

where r[ = y^x'^ + y'^ + z'^ , r'^ = ^J{x' - l)'^ + y'2 + ^'2 _ 



(13) 



df^ (1 + q)rl (1 + q)rl ' 
where ri = -^f x"^ -\- y'^ -\- z'^ , r2 — ^/(x — 1)^ + y'^ + z'^ . 

The equations of motion in the {x'S2y'z') coordinate system become: 

d^ d]/_^ ^' Q'jx' - 1) , 

dt^ dt 1 + q' (l + g'K^ {l + q>)r'^3 



df" dt ^ (i + g'K^ (i + g')r;3 ^ ^ 



(16) 



It can be easily verified that the equations of coordinate transformation are: 

x' = 1 — x , y' = y 1 z = z . (17) 

5. 'Similctr' initial conditions 



We have three differential equations of motion of second degree, therefore the initial 
conditions consist in three initial positions and three initial velocities. But in the following 
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we shall consider the motion of the test particle in the orbital plane, so the equations of 
motion will be (11)-(12) and (14)-(15) respectively. By consequence we need two initial 
positions and two initial velocities. 

5.1. 'Similcir' initial positions 

Considering Pq the point corresponding to the initial position, we denote Pq the 
projection of this point to the abscissa's axis. 

Following the same idea as in section 3, we define the initial abscissa as a given 
number whose absolute value represents the distance of Pq to the more massive star, 
measured in a given sense of abscissa 's axis. 

We define the initial ordinate as a given number representing the distance PqPo (see 
Figure 3, where Si was considered the most massive star). 

By consequence we have two 'similar' initial positions: Po(xo,yo) a-nd PQ{xQ,y'Q), where 

xq^i + xo , y'o = yo ■ (18) 

5.2. 'Similcir' initial velocities 

We define the initial velocity using the following three conditions: 

1. The components of initial velocity have given absolute value. 

2. The abscissa's axis, the initial velocity vector, and the Keplerian angular velocity 
vector form a trihedron with a given orientation (positive or negative). 

3. The angle formed by the initial velocity vector and the positive abscissa 's semi-axis 
has a given type (acute or obtuse). 



Then, for Vo{vox,voy) and VQ{vQ^,VQy) we have: 

'^'OX = -^OX , VQy = -VQy 



(19) 



Remark: All three conditions are necessary. If only (i) is considered, we have four 
possibilities for Vq (see Figure 4 (a)). If the conditions (i) and (ii) are considered, there are 
two possibilities for Vq (see Figure 4 (b), where a positive orientation of the trihedron is 
taken). If all the three conditions are considered, there is only one possibihty (see Figure 
4 (c) for an acute angle formed by the initial velocity vector and the positive abscissa's 
semi-axis) . 



6. 'Similar' trajectories 

As a numerical application in Figure 5, there are given two 'similar' trajectories 
for the binary system Earth-Moon { q — 0.0121 and q' — 82.45). The initial 
conditions are: Xq = —0.2 , uq = —0.5 , Vq^ = —0.8 , VQy = —0.6 , respectively 
Xq — 0.8 , Uq — —0.5 , Vq^ — 0.8 , — 0.6 . The time of integration is one 
Keplerian period. In Figure 6, the 'similar' trajectories for the same binary system, 
with initial conditions: xq = 0.6 , yo = 0.4 , vq^ = 0.5 , v^y = , respectively 
Xq = 1.6 , 7/q = 0.4 , v'q^ = —0.5 , VQy = are presented. The time of integration is 
one Keplerian period. In Figure 7, there are given the 'similar' trajectories for the same 
conditions as in Figure 6, but the time of integration is 200 Keplerian periods. 
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Fig. 4. — 'Similar' initial velocity 
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Fig. 5. — 'Similar' trajectories . 




Fig. 6. — Others 'similar' trajectories 
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7. The use of 'similar' coordinate systems to demonstrate some analytical 

relations in the Roche geometry 

We shall show now the consequence of 'similar' coordinate system use on the 
equipotential curves (see Figure 8), and on the equilibrium points' positions, in the frame 
of the Roche geometry. 

From equations of motion (11)-(12) we obtain the equation of potential function: 

2 



X 



+ y + 



+ 



2q 



l + qj ' ' (l + g)ri ' (l + g)r2 
(We prefer to use this form, because in the sec o nd p art of this paragraph we shall 



demonstrate some analytical relations of 



Seidovl ( I2OO4J ). who denoted the expression in 



square brackets with ip{x,y), (see equation 21).) Similarly, from equations (14)-(15), the 
corresponding potential function is: 

2 o 



'^'{x',y' 



1 



X 



1 



l + q' 



2q' 



1 + q')r[ (1 + q')r'^ 



In Figure 8 are represented the equipotential curves for a binary system characterized 
by the mass ratio q = 0.2 (and by consequence q' = 5). Let us consider Lj(xj,0), i = 1,5 
the Lagrangian points. For the first La grangian point w e obtained xi = 0.65856 and 



3j 1 



0.34144. The Jacobian constant ( (jSzebehely 



19671) p. 16 is Ci = C[ = 3.74897. For 



the second Lagrangian point, X2 = 1.43808 and X2 = —0.43808 and C2 = C2 = 3.53634. 
For the third Lagrangian point, X3 = —0.90250 and X3 = 1.90250 and C3 = C3 = 3.16504. 
The points L4 and L5 form equilateral triangles with Si and 6*2. The numerical results 
obtained above are normal, because x' = 1 — x (see equation(17)). In what concerne 



a = c' 



1, 3, th e equa. 



Jacobian constant. In 



i ties a re normal if we think at the physical meaning of the 



Seidovl (I2OO4I ). there are obtained some analytical relations in the 



Roche geometry and one of them can be considered as an analytical demonstration of the 
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relation Ci = C'^ , « = 1, 3. He obtained a very important correlation between the potential 

and the mass ratio (g) on one side, and the Lagrangian points Li(xi,0) and L2(x2,0) 
on the other side, in the frame of the classical Roche model. The relations obtained by 
Seidov are: 

^pi{q,xi) = ipi (1/g, 1 - Xi) , 

q{xi) ■ q{l - Xi) = I , i = T;2. (20) 

The first relation of (20) correspond to Cj = , i = 1,2 . 

This paragraph will show an elegant and unitary proof for formulae (20), using the 
'similar' coordinate systems. Using q, q', ri, r[, r^, r'^ we will obtain these relations for 
i = 1,3 and we will expand the expressions concerning the potential for i = 4, 5. So, the 
analytical relations of Seidov will be analyzed for all Lagrangian points, in a very easy 
manner. 



The equation of the potential given in 



Seidovl ( I2OO4J ) in the {xS\y) system is: 



^{x,y) = [x-j^] +y^ + 



2 2q 

(21) 



(l + g)ri (l + g)r2 

and if we use 'similar' coordinate systems, the equation of the potential in {x'S2y') system 
becomes: 



^l;\x',y') = [x'-^] +y 



where 



ri 



+ ^ + . (22) 

{l + q')r[ (l + g0^2 



a/x^ + ?/2 ^ r2 = a/(1 - X)2 + y2 



y'\ 
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Li{xi,0) being equilibrium points, we have dip/dxi = and 



dip' /dx[ = 0, i = 1, 3 (see Figure 9). From these relations we obtain: 

(1 -xif{l + xi + xl) 



q{xi) 
q{x2) -- 



xf (3 — 3xi + x\) 

{l-x\f{l + x\+x'^) 
xf(3-3x'i + x;2) ' 

(X2 - 1)^(1 +X2 
Xl{2 - X2)(l - X2+ Xfj 

(l-4)^(xf + 1) 
(l-X3)^(xi + l) 

x|(3 - 3x3 + xfj 
(l-4)2(xf-l) 



Xf(2-X'3)(l-X'3 + X?) • 

From (17) one can observe that = 1 — Xj , i = 1,3, and therefore g(xj) = ^7(^3^; 



(23) 

(24) 
(25) 
(26) 
(27) 
(28) 



and q(x,; ) ■ q'jl — Xj) 



given m 



1 ,2 = 1,3. These relations can be compared with the relations (20) 



Seidovl mm . for i = 1,2. 



For i = 4, 5 we have ^), L^{\, -^), ((Royl 



1988 



Murray! l2005l ) . (the coordinates 



of L4 and L5 are independent of g), therefore q is not a function of Xj, z = 4, 5. Here 
ri = 1 = r2. 

From (20) and (21), the potential corresponding to the Lagrangian points are: 

2 



tp{q,xi) 



Xi 



l + q 



+ 



+ 



2q 



(l + g)xi (l + g)(l-xi) 



(29) 



ilj\q,x\] 



X 1 



1 + g' 



+ 



y + 



2g' 



(l + g')a;'i (i + g')(i -2;'i) 



(30) 
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X2 



l + q 



+ 



2q 



il + q)x2 (l + g)(l-X2) 



(31) 



l + q' 



(1 + q')x'2 (1 + q')il - X 



2q' 



(32) 



X3 



l + g 



+ 



2q 



(l + g)x3 (l + g)(l-X3) 



(33) 



'^'i.q^x'^) 



2q' 



{l + q')x'^ {l + q'){l-x'^) 
Using the equations: x'j = 1 — Xj, and q' = 1/q we obtain: 



i>{q,Xi) = i>'{l/q, 1 - X,) 



i = 1,3 . 



These relations can be compared with relations (20) given in ISeidovi ( 12004J ) for i = 1,2. 



(34) 



For -^4(1, ^), using the equations (20) and (21) we obtain 



3 + 5g + 3g2 



ip'iq') 



3 + 5g' + 3q 



12 



(l + g)2 ' (l + g')' 

and because q' = 1/q, for L4 we obtain ?/;(g) = '?/''(l/g). 

For L5(|, — ^) the same equation is obtained. 

From equations (23)- (28) and (29)- (34) we obtain analytical formulae for the potential 
as a function of Lagrangian point positions: 

Ax\ - IQxl + lAx\ + 14xf - 41x^ + 40xf - 21x1 + 12xi - 3 



{x\ - 2x1 + 2xi - ly 
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Axl - Uxl + 18x1 + 94 - 36x| + 27x1 - 4x2 - 1 



- 2x1 + 2a:2 - 1)^ 



Axl 



lAxl + 18x^ 



29xf + 40xf 



27x1 + 12x3 - 3 



2xi + xl- 2x3 + 1)^ 



16x; ^ + 14x; 6 + I4x; 5 - 41x; ^ + 40x; ^ _ 21x[ 2 + 12x; - 3 



X 



/ 4 



2x; 



/ 3 



X' 



/ 2 



+ 2x; - 1) 



Ax' ^ - 14x', 6 + 18x; 5 - 29x; ^ + 40x; 3 - 27x; 2 + I2x' 



x: 



/ 4 



2x'2 3 + x'2 2 - 2x'2 + 1)^ 



4x'3 7 - 14x'3 6 + 18x'3 5 + 9x3 



/ 4 



36x'3 + 27x'3 2 - 4x'3 



X 



/ 4 



2x3 ^ + x; 



'3 ' + 24 



The equation for 'ip{x2) is obtained also by Seidov (see equation (9) in 



Seidovl mm '] 



8. Conclusion 

The 'similarity' relation defined in section 2 belongs to equivalence relations' family. 
For the time being it has only a theoretical value, completing classical method of study 
of the restricted three-body problem (see sections 4, 5, 6), and allowing for a more 
elegant demonstration of some analytical relations from the geometry of the Roche model 
(section 7). The use of 'similar' coordinate systems imposes the typical definitions of mass 
ratio, of distances from the test particle to the components of the binary system, and of 
initial conditions necessary to integrate the differential equations of motion. The 'similar' 
trajectories are not like- wise, but have a similar topology. 



The use of 'similar' coordinate s ystems 



the analytical relations obtained by 



helped us to create an elegant and easy proof for 



Seidovl ( I2OO4I ). To close the circle, we have completed 



these relations by analyzing the problem of mass ratio and potential, as function of 
Lagrangian point positions for all five Lagrangian points. So, the use of 'similar' coordinate 
systems in the restricted three-body problem enable us to complete the study of the Roche 
geometry with some new elements. 
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Fig. 8. — The 'similar' equipotential curves for the Roche model 
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Fig. 9. — The 'similar' coordinate systems and the Lagrangian points Li, L2, and L3. 
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